clear all; close all; clc;

R=6400000;
h=19100000;
Ror=R+h;
Tobr=12*3600;
dT=Tobr/1000;
t=(0:dT:Tobr);
phi_pred=tan(R/Ror)*180/pi;%предельное значение по углу
phi=(t/Tobr)*2*pi;%изменение угла
dphi=(dT/Tobr)*2*pi*180/pi;
N=phi_pred/dphi;
V1=2*pi*Ror/Tobr;

x = Ror * sin(phi);
y = Ror * cos(phi);

figure(1);
plot(x, y)

vx = diff(x)/dT;
vy = diff(y)/dT;
figure(2)
v2=(vx.^2+vy.^2).^(1/2);
plot(t(1:1000), v2)


xr=x;
yr=y-R;

Vx2=xr(1:1000).*vx;
Vy2=yr(1:1000).*vy;
V2=(Vx2.^2+Vy2.^2).^(1/2);

rv=(xr.^2+yr.^2).^(1/2);

plot(phi(1:1000)*180/pi,V2(1:1000)./rv(1:1000))

